Temporal disorder in up-down symmetric systems 



(N 

o 

(N 

P-h 

(N 



o 

S 

c3 



O 

o 



> 
m 
O 
in 
in 

o 



X 



Ricardo Martfncz-Garci'a, 1 Fcderico Vazquez, 2 Cristobal Lopez, 1 and Miguel A. Munoz 3 

1 IFISC, Instituto de Fisica Interdisciplinar y Sistemas Complejos (CSIC-UIB), E-07122 Palrna de Mallorca, Spain 
2 Max- Planck- Institut fur Physik Komplexer Systeme, Nothnitzer Str. 38, 01187 Dresden, Germany 
3 Instituto Carlos I de Fisica Teorica y Computacional, Universidad de Granada, 18071 Granada, Spain 

(Dated: February 27, 2012) 

The effect of temporal disorder on systems with up-down Z2 symmetry is studied. In particular, we 
analyze two well-known families of phase transitions: the Ising and the generalized voter universality 
classes, and scrutinize the consequences of placing them under fluctuating global conditions. We 
observe that variability of the control parameter induces in both classes "Temporal Griffiths Phases" 
(TGP). These recently-uncovered phases are analogous to standard Griffiths Phases appearing in 
systems with quenched spatial disorder, but where the roles of space and time are exchanged. TGPs 
are characterized by broad regions in parameter space in which (i) mean first-passage times scale 
algebraically with system size, and (ii) the system response (e.g. susceptibility) diverges. Our results 
confirm that TGPs are quite robust and ubiquitous in the presence of temporal disorder. Possible 
applications of our results to examples in ecology are discussed. 

PACS numbers: 64.60.Ht, 05.40.Ca, 05.50.+q, 05.70.Jk 



I. INTRODUCTION 



Systems with up-down Zi symmetry -such as the Ising 
model- are paradigmatic in Statistical Mechanics. Some 
of them -such as the voter model- exhibit absorbing 
states, a distinctive feature of non-equilibrium dynamics 
[lHH . Absorbing states are configurations of the system 
characterized by the lack of fluctuations, where the dy- 
namics becomes frozen and the system remains trapped. 
In the last years, a great interest has been given to this 
class of models with two symmetric absorbing states [J- 
Q , which are of high relevance in diverse problems in the 
ecological, biological, and social sciences, such as species 
competition neutral theories of biodiversity .al- 
lele frequency in genetics opinion formation m 1 
epidemics propagation 



14( , or language spreading 



Phase transitions into absorbing states are quite uni- 
versal. Systems exhibiting one absorbing state belong 
generically to the very robust Directed Percolation (DP) 
universality class and share the same set of critical expo- 
nents, amplitude ratios, and scaling functions. When this 
general rule is broken it is because some extra symmetry 
or conservation law is present in the system [JQ. This 
is the case of the class of systems with ^-symmetric 
absorbing states, which may exhibit a phase transition 
with critical scaling differing from DP, usually referred 
as Generalized voter (GV), also called "parity conserv- 
ing", "DP2" or "directed Ising", universality class (see 
PiS Q and references therein) . More interestingly, ana- 
lytical and numerical studies 0-0 have shown that, 
depending on some details, .^-symmetric models may 
undergo either a unique GV-like phase transition separat- 
ing an active/symmetric phase from an absorbing one or, 
alternatively, such a transition can split into two separate 
ones: a Ising- like transition in which the Z^ symmetry is 
broken, and a second DP-likc transition below which the 
broken-symmetry phase collapses into the corresponding 
absorbing state. In particular, a general Langcvin equa- 



tion, aimed at capturing the phenomenology of these sys- 
tems, was proposed in [4|; depending on general features 
it may exhibit a DP, an Ising, or a GV transition. 

In many situations, Zi symmetric systems are not iso- 
lated but, instead, affected by external conditions or by 
environmental fluctuations. The question of how exter- 
nal variability affects diversity, robustness, and evolution 
of complex systems, is of outmost relevance in different 
contexts. Take, for instance, the example of the neutral 
theory of biodiversity: if there are two ^-symmetric (or 
neutral) species competing, what happens if depending 
on environmental conditions one of the two species is fa- 
vored at each time step in a symmetric way? Does such 
environmental variability enhance species coexistence or 
does it hinder it? [H-io|. 

Motivated by these questions, we study how basic 
properties of Z^ symmetric systems, such as response 
functions and hrst-passage times, are affected by the 
presence of temporal disorder. 

Some previous works have explored the effects of fluc- 
tuating global conditions in simple models exhibiting 
phase transitions [2lM23l ]. Temporal disorder has been 
shown to be a highly relevant perturbation around DP 
phase transitions in all dimensions (in apparent contra- 
diction with the Harris criterion for the relevance of dis- 
order [U|), while temporal disorder has been shown to 
be relevant at Ising transition only at and above three di- 
mensions. More recently, a modified version of the sim- 
plest representative of the DP class -i.e. the Contact- 
Process- equipped with temporal disorder was studied 
in [24j]. In this model, the control parameter (birth prob- 
ability) was taken to be a random variable, varying at 
each time unit. As the control parameter is allowed to 
take values above and below the transition point of the 
pure contact process, the system alternates between the 
tendencies to be active or absorbing. As shown in [24| 
this dynamical frustration induces a logarithmic type 
of finite-size scaling at the transition point and gencr- 
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ates a subregion in the active phase characterized by a 
generic algebraic scaling (rather than the usual exponen- 
tial, Kramers-like, behavior) of the extinction times with 
system size. More strikingly, this subregion is also char- 
acterized by generic divergences in the system suscepti- 
bility, a property which is reserved for critical points in 
pure systems. This phenomenology is akin to the one 
in systems with quenched "spatial" disorder [26| . which 
show algebraic relaxation of the order parameter, and sin- 
gularities in thermodynamic potentials in broad regions 
of parameter space: the so-called, Griffiths Phases [25| . 
The remarkable peculiarities of standard Griffiths phases 
stem from the existence of (exponentially) rare -locally 
ordered- regions which take a (exponentially) long time 
to decay, inducing an anomalously slow decay in the dis- 
ordered phase. 

In the case of temporal disorder, an analogy with Grif- 
fiths phases can be made, in the sense that very long 
(exponentially rare) time intervals (corresponding to an 
absorbing phase of the pure model) of the control pa- 
rameter have a large influence on the system dynamics 
even when the overall system is in its active phase. These 
phenomenological similarities between systems with spa- 
tial and temporal disorder led to introduce the concept 
of "Temporal Griffiths Phases" (TGP) 0. 

In order to investigate whether the anomalous behavior 
that leads to TGPs around absorbing state (DP) phase 
transitions is a universal property of systems in other 
universality classes -and in particular, in up-down sym- 
metric systems- we study the possibility of having TGPs 
around Ising and GV transitions. We scrutinize simple 
models in these two classes and assume that the cor- 
responding control parameter changes randomly in time, 
fluctuating around the transition point of the correspond- 
ing pure model, and study the susceptibility as well as 
mean-first passage times. We mainly focus on the mean- 
field (high dimensional) limit, since it allows for analyti- 
cal treatment via a Langevin approach, but we also pro- 
vide numerical results and some theoretical considera- 
tions for low dimensional systems. 

The paper is organized as follows. In section [TTl we 
develop a general mean-field description of models with 
varying control parameters in terms of collective vari- 
ables. In section Hill and HVl we show analytical and nu- 
merical results for the Ising and GV transitions, respec- 
tively. In section [Vj a short summary and conclusions 
are presented. 



II. MEAN-FIELD THEORY OF Z2-SYMMETRIC 
MODELS WITH TEMPORAL DISORDER 

Interacting particle models evolve stochastically over 
time. A useful technique to study such systems is the 
mean-field (MF) approach, which implicitly assumes a 
"well-mixed" situation, where each particle can interact 
with any other, providing a sound approximation in high 
dimensional systems. One way in which the mean-field 



limit can be seen at work is by analyzing a fully con- 
nected network (FCN), where each node (particle) is di- 
rectly connected to any other else, mimicking an infinite 
dimensional system. 

In the models we study here, states can be labeled 
with occupation-number variables pi taking a value 1 if 
node i is occupied or if it is empty, or alternatively by 
spin variables Si = 2pi — 1, with Si = ±1. Using these 
latter, the natural order parameter is the magnetization 
per spin, defined as 



m 
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where N is the total number of particles in the system. 
The master equation for the probability P(m, t) of having 
magnetization to at a given time t, is 

P(m,t + 1/N) = u) + (m-2/N,b)P(m-2/N,t) (2) 
+ w_ (to + 2/N, b)P(m + 2 fN, t) 
+ [1 — u)-(m, 6) — cj_|_(m, 6)] P(m, i), 

where w± (771,6) are the transition probabilities from a 
state with magnetization to to a state with magneti- 
zation to ± 2/N. This describes a process in which a 
"spin" is randomly selected at every time-step (of length 
dt = 1 /N) , and inverted with a probability that depends 
on m and the control parameter 6. The allowed magneti- 
zation changes in an individual update, A?n = ±2/N, are 
infinitesimally small in the N — > 00 limit. In this limit, 
one can perform a standard Kramers-Moyal expansion 
[29l [3fjj leading to the Fokker-Planck equation 

8P(m t) 8 1 8 2 

— m =~d^ [/(m ' 6)p(m ' t)]+ 2 [5(m ' &)p(m ' t)] ' 

(3) 

with drift and diffusion terms given, respectively, by 

/(m,6) - 2[w + (m J 6)-w_(m,6)], (4) 



5(777,6) 



4 [u+ (to, b) + uj- (to, 6)] 

N ' 



(5) 



From Eq. ([3]), and working in the Ito scheme (as justified 
by the fact that it comes from a discrete in time equation 
[2q|). its equivalent Langevin equation is (30| 



777 = /(m, 6) + y/g{m,b)r){t), 



(6) 



where the dot stands for time derivative, and 77 (t) is 
a Gaussian white noise of zero-mean and correlations 
(r](f)r)(t')) = S(t — t'). The diffusion term is propor- 
tional to 1/y/N, and therefore, it vanishes in the ther- 
modynamic limit (N — > 00), leading to a deterministic 
equation for to. 

The drift and diffusion coefficients in Eq. depend 
not only on the magnetization, but also on the parameter 
b. To analyze the behavior of the system when b changes 
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FIG. 1: (a) Typical realization of the colored noise £(t), a 
step like function that takes values between +1 and — 1. (b) 
Stochastic control parameter b(t) = bo + cr£(t) according to 
the values of the noise in (a), bo = 1-03 and a = 0.4. 



t/3. Then, we combine the two white noises into an 
effective Gaussian white noise, whose square amplitude 
is the sum of the squared amplitudes of both noises [3(| , 
and finally arrive at 



th = /o(m) + yjg (m) + Kj 2 (m) j(t), (10) 

where ( 7 (f)) = and (7(t)7(*')> = S(t - t'). 

In the next two sections, we analyze the dynamics 
of the kinetic Ising model with Glauber dynamics and 
a variation of the voter model (the, so-called, q-voter 
model) -which are representative of the Ising and GV 
transitions respectively- in the presence of external noise. 
For that we follow the strategy developed in this section 
to derive mean-field Langevin equations and present also 
results of numerical simulations (for both finite and infi- 
nite dimensional systems), as well as analytical calcula- 
tions. 



randomly over time, and following previous works [19i . 
|24| . we allow b to take a new random value, extracted 
from a uniform distribution, in the interval (6o — °": bo+a) 
at each MC step, i.e., every time interval r = 1. Thus, 
we assume that the dynamics of b(t) obeys an Ornstein- 
Uhlenbeck process 



b(t) = b + a£(t), 



(7) 



where is a step- like function that randomly fluctuates 
between —1 and 1, as depicted in Fig. QJi. Its average 
correlation is 

|(1 - |At|/r) for |At| < t 
<£(i)£(f + At)>=< (8) 
for \At\ > t, 

where the bar stands for time averaging. Parameters &o 
and a are chosen with the requirement that b takes val- 
ues at both sides of the transition point of the pure model 
(see Fig. \Tjp), that is, the model with constant b. Thus, 
the system randomly shifts between the tendencies to be 
in one phase or the other (see Fig. The model presents 
both intrinsic and extrinsic fluctuations, as represented 
by the white noise n(t) and the colored noise re- 
spectively. Plugging the expression Eq. for b(t) into 
Eq. (HU), and retaining only linear terms in the noise one 
readily obtains 



m = f (m) + \Jga(m) r)(t) + j (m) £(t), 



(9) 



where fo(m) = f(m,b ), g {m) = g(m 7 b ) and jo(m) is 
a function determined by the functional form of f(m, b), 
that might also depend on bg. To simplify the analysis, 
we assume that relaxation times are much longer than the 
autocorrelation time r, and thus take the limit r — > 
in the correlation function Eq. ([5]), and transform the 
external colored noise £ into a Gaussian white noise with 
effective amplitude K = < £(£)£(* + At) > dAt = 



III. ISING TRANSITION WITH TEMPORAL 
DISORDER 

We consider the kinetic Ising model with Glauber dy- 
namics [HI, as defined by the following transition rates 



1 — Si tanh 



2d 



(11) 



The sum extends over the 2d nearest neighbors of a given 
spin i on a d-dimcnsional hypercubic lattice, and b = 
J/3 is the control parameter. J is the coupling constant 
between spins, which we set to 1 from now on, and /3 = 
(/cbT) -1 . Note that b in this case is proportional to the 
inverse temperature. 



A. The Langevin equation 

In the mean-field case, the cubic lattice is replaced by 
a fully-connected network in which the number of neigh- 
bors 2d of a given site is simply N — 1. Then, the tran- 
sition rates of Eq. (| 1 1 [) can be expressed as 

n ± (m,b) = n(=F ->-±) = - [1± tanh (6m)] . (12) 

which implies u±(m, b) = 1= F m tt±(m, 6) for jumps in the 
magnetization. Following the steps in the previous sec- 
tion, and expanding Q± to third order in m, we obtain 



m = aom—com 



1 — bom 2 
N 



Ka 2 m 2 (l-b 2 m 2 ) 2 y(t), 
(13) 



where bo is the mean value of the stochastic control pa- 
rameter, a = bo — 1, and c = 6q/3. 

ao , 

2'" 1 4 

the deterministic term of Eq. (fTB")) has the standard shape 



The potential V(m) = — m? + ^-m 4 associated with 
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FIG. 2: Schematic representation of the fluctuating control 
parameter in the Ising model with Glauber dynamics. It 
makes the system shift between the ordered phase to the dis- 
ordered one. 



of the Ising class, that is, of systems exhibiting a sponta- 
neous breaking of the Z 2 symmetry. A single minimum 
at m = exists in the disordered phase, while two sym- 
metric ones, at a/c exist below the critical point. 



B. Numerical Results 

By performing simulations of the Langevin equation 
Eq. ([15]) or, equivalently, by performing Monte Carlo 
simulations of the particle system on a fully connected 
network, we observe that the critical point is located at 
the value &o,e = 1 (identical to that of the pure case, 
b c ,pure = 1; see below for details). Having located the 
transition point, we study two magnitudes that are rel- 
evant in systems with temporal disorder (24|, the mean 
crossing time (or mean-first passage time) and the sus- 
ceptibility. The crossing time is the time employed by the 
system to reach the disordered zero-magnetization state 
for the first time, starting from a fully ordered state with 
|m| = 1 (sec Fig. [3]). Crossing times were calculated by 
numerically integrating Eq. (|13[) for different realizations 
of the noise 7 and averaging over many independent re- 
alizations. Results are shown in the Fig. 2J 

As it is characteristic of TGPs 
point, the scaling is of the form T ' 



24 



at the critical 
[\nN] a with a nu- 
merically determined exponent value a = 2.81 for a = 0.4 
(see inset of Fig.[4]). This is to be compared with the stan- 
dard power-law scaling T ~ N 13 characteristic of pure 
systems, i.e. for a = 0. Moreover, a broad region show- 
ing algebraic scaling T ~ N s with a continuously varying 
exponent 5 (bo) (8 — > as bo — > 6j c ) appears in the or- 
dered phase bo > bo_ c . Both a and 5 are not universal 
and depend on the noise strength a. Finally, in the disor- 
dered phase the scaling of T is observed to be logarithmic, 
T~lniV. 

We have also performed Monte Carlo simulations of 
the time-disordered Glauber model on two- and three- 
dimensional cubic lattices with nearest neighbor inter- 
actions. In d — 2, a shift in the critical point was 
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FIG. 3: Single realization of the stochastic process. The sys- 
tem starts with all the spins in the same state (m = 1) and 
the dynamics is stopped when it crosses m = 0, which de- 
fines the crossing time in the Ising model. We take a — 0.4. 
bo — 0.98 and system size N = 10 6 . On the right margin we 
sketch the mapping of the problem to a Random Walk with 
jump length \Am\ = 2/N. 
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FIG. 4: Main: Log-log plot of the crossing time T(N) for 
the Ising Model with Glauber dynamics in mean field. Monte 
Carlo simulations on a FCN (circles) and numerical integra- 
tion of the Langevin equation Eq. (|13[) with a — 0.4 (solid 
lines). There is a region with generic algebraic scaling of 
T(N) and continuously varying exponents. Inset: log-log plot 
of T(N) vs. In AT. It is estimated at criticality T ~ (IniV) 2 - 81 . 



found: from b CiVure = 0.441(1) in the pure model to 
60, c = 0.605(1) for a = 0.4. However, the scaling be- 
havior of T with N resembles that of the pure model, 
with T ~ N@ at criticality (with an exponent numerically 
close to that of the pure model Q), and an exponential 
growth T ~ cxp(a/V), where c is a positive constant, in 
the ordered phase (Arrhenius law). Thus, no region of 
generic algebraic scaling appears in this low-dimensional 
system. On the contrary, in d = 3, results qualitatively 
similar to mean- field ones are recovered (see Fig. [5]). The 
critical point is shifted from b CtPure = 0.222(1) (calcu- 
lated in [3l|) to bo.c = 0.413(2), with a critical exponent 
a(d = 3) = 5.29 for a = 0.4, and generic algebraic scal- 
ing in the ordered phase. In conclusion, our numerical 
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FIG. 5: Main: Log-log plot of the crossing time T(N) for the 
Ising Model with Glauber dynamics in d — 3. Monte Carlo 
simulations on a regular cubic lattice with a = 0.4. There is a 
region b 6 [0.42, 0.46] with generic algebraic scaling of T(N) 
and continuously varying exponents. Inset: log-log plot of 
T(N) vs. ln(JV). It is estimated at criticality T ~ (IniV) 5 ' 29 . 



FIG. 6: Main: Log-log plot of the susceptibility as a func- 
tion of the external field for different values of bo 6 [6o,c — 
(j 2 /2, &o,c + o" 2 /2] with a = 0.1, obtained by integrating 
Eq. (|15jl . Generic divergences with symmetric exponents 
around the critical value bo. c = 1 are observed. 



studies suggest that the lower critical dimension for the 
TGPs in the Ising transition is d c = 3. This is in agree- 
ment with the analytical finding in p2j , establishing that 
temporal disorder is irrelevant in Ising-like systems be- 
low three dimensions. This result is to be compared with 
d c = 2 numerically reported for the existence of TGPs 
in DP-like transitions [24[ (observe, however, that tem- 
poral disorder, in this case, affects the value of critical 
exponents at criticality in all spatial dimensions). Fur- 
ther studies are needed to clarify the relation between 
disorder-relevance at criticality and the existence or not 
of TGPs. 

We have also measured the susceptibility Xi defined 
as the response function to an external field h in the 
vanishing field limit 



X 



dm 
fc->o oh 



(14) 



where fh denotes the stationary magnetization averaged 
over many independent realizations. In the presence of 
an external field, the transition rates become il±(m, b) = 
2 [1 ± tanh (b m + h)] . Expanding the hyperbolic tan- 
gent up to third order in m and to first order in h, we 
obtain the following Langcvin equation 



m = aom — com 3 + h (l — 6 2 m 2 ) +vKam(l — b m ) 7(t), 

(15) 

where we have considered the N — > oo limit (g = 0). 

The average magnetization fh for a given field h was 
calculated by integrating the Langcvin equation and then 
taking averages over noise realizations. The susceptibil- 
ity can be computed, for different values of bo, as the 
derivative of fh with respect to h. Generic divergences of 
the form \ ~ /i^+Constant (with v < as h — > 0) appear 
in a broad region bo <G [6o,c — c 2 /2, b 0:C + c 2 /2], centered 
around &o,o with symmetric exponents around the criti- 
cal point (see Fig. [5]). These results agree with those ob- 
tained through Monte Carlo simulations on a FCN (not 



shown). In finite dimensions, given the required large sys- 
tems sizes and small fields, we could not conclude about 
the existence or not of generic divergences. 



C. Analytical results 

Let us consider the Langevin equation Eq. (|13[) in the 
Thermodynamic limit (go( m ) = 0)- Given that the re- 
maining intrinsic noise comes from a transformation of a 
colored noise into a white noise, the Stratonovich inter- 
pretation is to be used to obtain its associated Fokkcr- 
Planck equation (see e.g. [28|) 



dP(m,t) 
dt 




K 

fo(m) + yio(m)i (m) 



P(m,t) 



2 dm 2 



[Kjl{m)P{m,t)). 



(16) 



Imposing the detailed balance (fluxless) condition, it is 
straightforward to obtain the steady state solution 



P s t{m) oc exp 



Com 



|2a /<7 -1 



(17) 



with a power-law singularity at the origin, which is a dis- 
tinctive trait of a Langevin equation with linear multi- 
plicative noise 1321 13311 . By performing a calculation anal- 
ogous to that in [24|, we have analytically computed the 
system susceptibility and found that x ~ h v + Constant, 
as mentioned earlier, and in agreement with previous re- 
sults found in [H HI [H, with 



(18) 



This, in particular, implies that the susceptibility di- 
verges when v < as h — > or, in terms of the control 
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parameter bo = 1+ao, when bo takes a value in the region 
1 — cr 2 /2 < bo < 1 + er 2 /2 centered at the critical point 
bo.c = 1- The values of the exponent v agree well with 
those of Fig. IS] at some distance from the critical point. 
For instance, an analytical value v = —0.40 for bo = 1.003 
corresponds to a numerical value v nU m = —0.39, and 
v = —0.60 for bo = —1.002 to a value v num = —0.59. 
However, the analytical exponent v = — 1 at the critical 
point is not in good agreement with the numerical result 
Vnum = ~0.88, indicating that the asymptotic regime has 
not been numerically reached. 

We next provide analytical results for the crossing 
time. Starting from the N-indcpendcnt Fokkcr-Planck 
equation Eq. (|16|) . an effective dependence on N is im- 
plemented by calculating the first-passage time to the 
state m = \2/N\ rather than to = 0. This is equivalent 
to the assumption that the system reaches the zero mag- 
netization state with an equal number N + = iV_ = N/2 
of up and down spins when |m| < 2/N. that is, when 
N/2 - 1 < N + < N/2 + 1. The mean- first passage time 
T associated with the Fokker-Planck equation Eq. (TTF 
obeys the differential equation [30j 



y f (m)T"(m) + 



K 

fo(m) + yjo(™)jo( m ) 



T'(m) 



= -1, 
(19) 

with absorbing and reflecting boundaries at \m\ = 2/N 
and \m\ = 1, respectively. The solution, starting at time 
t = from to = 1 is given by 



1 dy f 1 i>(z) 



T(m=l) = 2f * / 

J2/N TO Jv 



/N TO) Jy Kjo(z) 



(20) 



where 



[J 2/N K .lo( x ) 



Computing these integrals (see Appendix [A} we obtain 

(22) 



]n.N/(b - 1) for 6 < 1 
T~ I 3(lnN) 2 /a 2 for b = 1 



6(6 -l) 

N^~ for b > 1. 



These expressions qualitatively agree with the numerical 
results of Fig. @J showing that T grows logarithmically 
with N in the absorbing phase bo < 1, as a power law in 
the active phase bo > 1, and as a power of In AT (i.e. 
poly- logarithmically) at the transition point bo.c = 1- 
The critical exponent a — 2, as well as the exponents 
8 = 6(bo — l)/o" do not agree well with the numer- 
ically determined exponents. This is probably due to 
to the fact that we have neglected the 1/Vn term by 
taking go = 0, which becomes of the same magnitude 
as the jo term when \m\ approaches 2/N. Indeed, this 
was confirmed (not shown) by testing that analytical ex- 
pressions Eq. (|22[) agree very well with numerical inte- 
grations of Eq. (|13|) performed for go = 0, and setting 



the crossing point at m = 2/N. In summary, this an- 
alytical approach reproduces qualitatively -and in some 
cases quantitatively- the above reported non-trivial phe- 
nomenology. 



IV. GENERALIZED VOTER TRANSITION 
WITH TEMPORAL DISORDER 

We study in this section the GV transition which 
appears when a .^-symmetry system simultaneously 
breaks the symmetry and reaches one of the two absorb- 
ing states. A model presenting this type of transition is 
the nonlinear q- voter model, introduced in [35j ] . The mi- 
croscopic dynamics of this nonlinear version of the voter 
model consists in randomly picking a spin Si and flip- 
ping it with a probability that depends on the state of 
q randomly chosen neighbors of Si (with possible rep- 
etitions). If all neighbors are at the same state, then 
Si adopts it with probability 1 (which implies, in par- 
ticular, that the two completely ordered configurations 
are absorbing). Otherwise, Si flips with probability b. 
Therefore, the probability that Si flips at a given time, is 
a function of the fraction x of disagreeing (antiparallel) 
neighbors 



f(x, b) = x q + b[l - x q - (1 - x) q ] 



(23) 



where b is taken as the control parameter. Three types 
of transitions, Ising, DP and GV can be observed in this 
model depending on the value of q Hf|. Here, we focus 
on the q — 3 case, for which a unique GV transition at 
b c = 1/3 has been reported [H| . 



A. The Langevin equation 



In the MF limit (FCN) [H], the fractions of antiparallel 
neighbors of the two types of spins Si = 1 and Si = — 1 
are x = (1 — m)/2 and x = (l + m)/2, respectively. Thus, 
the transition probabilities are 



( u\ 1 T ™ / 1 ± TO 

w±{m,b) = — - — / I — - — ,b 



(24) 



Following the same steps as in the previous section, we 
obtain the Langevin equation 



m = — - — to(1 — to ) + 



(25) 



(l-?77 2 )(l + 6&o+m 2 ) 9K 



N 



+ — CT 2 TO 2 (1 - 7T1 2 ) 2 j(t). 



Let us remark that the potential in the nonlinear voter 
model (Fig. [TJ differs from that for the Ising model. Ow- 
ing to the fact that the coefficients of the linear and cubic 
term in the deterministic part of Eq. (|25[) coincide (ex- 
cept for their sign), the system exhibits a discontinuous 
jump at the transition point, where the potential min- 
imum changes directly from to = in the disordered 
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FIG. 7: Potential for the GV transition in a mean field ap- 
proach. The dashed, solid and dot-dashed lines correspond to 
the paramagnetic phase, critical point, and the ferromagnetic 
phase, respectively. 



phase to m = ±1 in the ordered one. Furthermore, the 
potential vanishes at the critical point 




FIG. 8: Main: Log-log plot of the ordering time as a function 
of the system size N in the MF q-voter model. Monte Carlo 
simulations on a FCN (dots) and numerical integration of the 
Langevin equation Eq. (|25p for values of b going from 0.330 
(bottom) to 0.370 (top), and a = 0.3 (lines). In the active 
phase a finite region with power law scaling is observed. Inset: 
log- log plot of T as a function of In N. At the critical point 
is T ~ [In AT] 3 ' 68 . 



B. Numerical Results 



C. Analytical results 



The ordering time, defined as the averaged time re- 
quired to reach a completely ordered configuration (ab- 
sorbing state) starting from a disordered configuration, 
is the equivalent of the crossing time above. We have 
measured the mean ordering time T by both, integrating 
the Langevin equation Eq. (|2"5j) and running Monte Carlo 
simulations of the microscopic dynamics on FCNs and fi- 
nite dimensions. In Fig. [8] we show the MF results. We 
observe that T has a similar behavior to the one found 
for the mean crossing time in the Ising model, and for 
the mean extinction time for the contact process [24| . 
That is, a critical scaling T ~ [In AT]" at the transition 
point &o.c = 1/3, with a critical exponent a = 3.68 for 
a = 0.3, a logarithmic scaling T ~ In N in the absorbing 
phase bo < bo lC , and a power law scaling T ~ N s with 
continuously varying exponent 6 (bo) in the active phase 
h> > b 0>c . 

Monte Carlo simulations on regular lattices of dimen- 
sions d = 2 and d = 3 revealed that there is no signifi- 
cant change in the scaling behavior respect to the pure 
model (not shown). The critical point shifts in d = 2 
and remains very close to its mean-field value in d = 3, 
but results are compatible with the usual critical (pure) 
voter scaling T^d ~ AT hxN and T^d ~ AT. In the absorb- 
ing phase T grows logarithmically with N , while in the 
active phase T grows exponentially fast with N , as in the 
pure-model case. Therefore, in these finite dimensional 
systems we do not find any TCP nor other anomalous ef- 
fects induced by temporal disorder, although we cannot 
numcricaly exclude their existence in d = 3. Such effects 
should be observable, only in higher dimensional systems 
(closer to the mean- field limit). 



The ordering time T can be estimated by assuming 
that the dynamics is described by the Langevin equa- 
tion Eq. (|25|) , and calculating the mean first-passage time 
from m = to any of the two barriers located at |m| = 1. 
It turns out useful to consider the density of up spins 
rather than the magnetization 



1 



(26) 



T is the mean first-passage time to p — starting from 
p = 1/2. The Langevin equation for p is obtained from 
Eq. (|25|) . by neglecting the 1/vN term and applying the 
ordinary transformation of variables (which is done em- 
ploying standard algebra, given that Eq. (|2"5*j) is inter- 
preted in the Stratonovich sense) is 



with 



P = A(p) + ^KC(p) 1 (t), 



A(p) = a p(2p-l)(l-p), 
C( P ) = 3<rp(2p-l)(l-p), 



(27) 



(28) 



where ao = 1 — 36o- 

Now, we can follow the same steps as in section IIII CI 
for the Ising model, and find the equation for the mean 
first-passage time T(p) by means of the Fokkcr-Planck 
equation. The solution is given by (see Appendix B) 



lnA/(36 - 1) for b < 1/3 
; (lnA073cr 2 for b = 1/3 

for b Q > 1/3. 



2(b -l/3) 

AT — ^ — 



(29) 
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These scalings, which qualitatively agree with the nu- 
merical results of Fig. [3] for the q-voter, show that the 
behavior of T is analogous to the one observed in the 
Ising transition of section IIIII and in the DP transition 
found in [24| . Therefore, we conclude that TGPs appear 



around GV transitions in the presence of external varying 
parameters in high dimensional systems. 

For the GV universality class the renormalization 
group fixed point is a non-perturbative one |37[, becom- 
ing relevant in a dimension between one and two. A field 
theoretical implementation of temporal disorder in this 
theory is still missing, hence, theoretical predictions and 
sound criteria for disorder relevance are not available. 



turcs. 



Acknowledgments 

R.M-G. is supported by the JAEPredoc program 
of CSIC. R.M-G. and C.L. acknowledge support from 
MICINN (Spain) and FEDER (EU) through Grant No. 
FIS2007- 60327 FISICOS. MAM acknowledges finan- 
cial support from the Spanish MICINN-FEDER under 
project FIS2009-08451 and from Junta de Andalucfa 
Proyecto de Excelencia P09FQM-4682. We are grate- 
ful to J. A. Bonachela for useful discussions and a critical 
reading of the manuscript. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the effect of temporal disorder 
on phase transitions exhibited by Zi symmetric systems: 
the (continuous) Ising and (discontinuous) GV transi- 
tions which appear in many different scenarios. We have 
explored whether temporal disorder induces Temporal 
Griffiths Phases as it was previously found in standard 
(DP) systems with one absorbing state. By performing 
mean-field analyses as well as extensive computer simu- 
lations (in both fully connected networks and in finite di- 
mensional lattices) we found that TGPs can exist around 
equilibrium (Ising) transitions (above d = 2) and around 
discontinuous (GV) non-equilibrium transitions (only in 
high-dimensional systems). Therefore, we confirm that 
TGPs may also appear in systems with two symmetric 
absorbing states, illustrating the generality of the under- 
lying mechanism: the appearance of a region, induced by 
temporal stochasticity of the control parameter, where 
first-passage times scale as power laws of the system size 
and where the susceptibility diverges. Temporal disorder, 
makes the ordered/active phase less stable and makes the 
system highly susceptible to perturbations. This appears 
to be a rather general and robust phenomenon. 

It also seems to be a general property that TGPs do 
not appear in low dimensional systems, where standard 
fluctuations dominate over temporal disorder. In all the 
cases studied so far, a critical dimension d c -at and below 
which TGPs do not appear- exist (d c = 1 for DP tran- 
sitions, d c = 2 for Ising like systems, and d c ~ 3 for GV 
ones). Calculating analytically such a critical dimension 
and comparing it with the standard critical dimension for 
the relevance/irrelevance of temporal disorder at the crit- 
ical point (i.e. at the renormalization group non-trivial 
fixed point of the corresponding field theory) remains an 
open and challenging task. 

A relevant application of our results is found in models 
of ecosystems. In this case, first-passage times are re- 
lated to typical extinction times, and studying how such 
extinction times are affected by system size (e.g. habitat 
fragmentation) is a problem of outmost relevance. Fu- 
ture research might be oriented to the effect of temporal 
disorder on the formation and dynamics of spatial struc- 



Appendix A: Analytical calculations of the crossing 
time for the mean field Ising model 

The mean first passage time to reach an absorbing bar- 
rier at |m| = 2/N starting from \m\ = 1 can be expressed 
as [30| . 



T(m = I) 



dV fl %*dz, (Al) 



/n ip(y) J y Kjl{z) 



with 



2/N 



which involves 6 th and 4 th order polynomial functions. 
In order to make the integral simpler, functions are ex- 
panded up to 3 rd order, 

fo(m) + yj ("i)io( m ) ~ m(r~sm 2 ) 

Kj$(m) « urn 2 , (A3) 

with uj = cr 2 /3, r = ao + oj/2, s = (cq + 2w6g). A second 
simplifying assumption is to take 1 as the lower integra- 
tion limit in Eq. (|A2[) instead of 2/N (justified because 
ip(z) appears both in the numerator and the denomina- 
tor of T(m) and the contribution of this limiting value is 
negligible). Therefore, Eq. (|A2[) becomes 

4>{z) = cxp f 2z ' {r ~* Z ' 2) dz' = z a e^- z2 \ (A4) 



where a = 2r/ui and f3 = s/uj. The first passage time is 
written as 



where it has been defined 

- 1 tf>(z) ^ rl 



(A5) 



[ Mrdz = °- I z°-^dz, (A6) 

Jy R Jo( z ) U Jy 



which presents a singularity when a = 1 (bo = 1 = 6o jC 
This case will be studied separately. 
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1. Case a / 1 



Integrating by parts Eq. (|A6|) . 



-P - e ~P% < y a ~ x 

a — 1 



2/3 



1 z a e-P z 
a- 1 



(A7) 

This integral can be solved again integrating by parts, 
and so on, recursively, 



1 ~ . i - c -/3(y 2 -i) y "-i+2fc 



nt a-l + 2i 



(A8) 



Therefore 
2 



CO 



(2/3) fe 



T = -E- , 



[-^(A) — h(k, A)] , (A9) 



where 



J X (A) 
/a(ft,JV) 



y-"^^ - 1] dy, 



l/N 
1 

1/JV 



(A10) 
(All) 



-Zi(A) is solved by parts. A recursive integration similar 
to the one in Eq. (|A6[) has to be performed, 



(-2/3)'[l- jV-*-i-2i e 0(i/" 2 -i)] 



z=o U J= o a- l + 2j 

On the other hand, l2(k,N) is easily solved 

„ N f -InfA- 1 ) = In A for fc = 0, 

The final expression for the first passage time is 



T = 2 fhW-hiN 



a — 1 

3 (2/3) fc [7 1 (A)-(l-7V- 2fc )/2fc] 



E 

fe=i 



n « =0 a-l + 2i 



(A14) 

whose asymptotic limit A — > oo has two different cases, 
a. a < 1 

In this case, a — 1 — 2/ < when I > so in I\ (A) 

1 _ N a-l-U f Ji(l/N*-l)„ lt (A15) 

which leads to 

(-2/3)' 



i= oii;-=oi+2j 



C(a,0). (A16) 



We have 
2 



r 



C(a,fl- In AT | (2/?) fc (C(a,/?)-(2fc)- 1 ) 

(A17) 

and finally, 



T 



uj(a — 1) 



b. a > 1. 



In A. 



(A18) 



Considering that A"" 1 > •A"" 1 " 2 ', V7 > 0, only the 
first term is relevant in Eq. (|A12[) for h(N). Then 



7i(A) 



l-c^A"" 1 e _/3 A a - 



1 - a 



1-Q ' 



(A19) 



and in the asymptotic behavior (A 3> 1) of the mean 
escape time 

T k, K(a, /3)A Q ~ 1 21nJV ^ jya-i_ ( A20 ) 

w(a — 1) 



2. Case a = 1. Critical point 



We need to solve 



< ai2 > '^'lM) d ' = ^l s " e ~" w (A21) 

Expanding the exponential function and integrating, it is 



(Ai3) !(y) 



■J 



-iny + E 



(_/?)*(! _2y) 2fe 



k=i 



k\2k 



(A22) 



Taking Eq. ([A"22]) into Eq. flX5 



2e^ 



where 



00 (— /T> fe 



k=l 



k\2k 



(A23) 



h{N) 
h(N) 
I 5 (k,N) 



Jl/N 



l/N 



l/N 
1 

l/N 



yk-l e P(v 2 -l) d y. 



(A24) 



First of all, let us consider the solution of /3(A) integrat- 
ing by parts, so that 



(In A) 2 



l/N 



(A25) 
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and wc obtain 



(In A) 2 



-V +2f3-0{N- 1 )+0 ( 



which scales in the asymptotic limit as 

, . (In TV) 2 a 
h{N) ~ ^-^-e^. 



In A 
(A26) 



(A27) 



Appendix B: Analytical calculations of the crossing 
time for the mean field nonlinear q-voter model 



After performing the change of variables of Eq. (|2(il) . 
the absorbing barrier is placed at p = 1/A and the re- 
flecting one at p = 1/2, (which is the initial point). The 
mean first passage time is given by 



On the other hand, the leading behavior when the size of 
the system is big enough (TV ^> 1) is 



h(N) ~ e-^lniV + C^). 



(A28) 



To solve the last integral, Is(k, A), the exponential func- 
tion has to be expanded as well. It is 



U(k + 2l) 



1 — TV 



-2l-k 



) ~ constant, 



(A29) 

when TV ^> 1. It finally leads to an expression for T at 
criticality 



T ! 



2e"^ 



! (lnTV) 2 



E 



til 

k\2k 



'\nN + C' 4 (P) 



(A30) 

In the limit of very large system sizes (TV 3> 1) the mean 
escape time scales as 



T 



(In A) 2 



1 OO 



k=l 



k\k 



]nN + K(/3), (A31) 



which asymptotically becomes 



(A32) 



Summing up, the time taken by the system for reaching 
m = 2/TV from an initial condition m = 1 is 



T 



uj(a-l) 
(In AT) 2 



A' 



In TV for a < 1, 

for a = 1, 
for a > 1. 



(A33) 



with 



i/jv Hv) Jy K 3o( z ) 



, 2A(z') + KC(z')C'(z') 
ip(z)=exp dz . (B2) 



KC 2 (z') 



l/N 



We expand the polynomials in Eq. (|B2[) up to second 
order, using Eq. (|28|) . it is 



A(p) + ^C0>)C"(p) 



p(r - sp), 

2 



(B3) 



where we have defined w = 3<r , r = ^ — 1 + 36o and 
s = 3r. These polynomials are similar to the ones ob- 
tained for the Ising model, but with redefined parame- 
ters. The integrals are done in a very similar way, and one 
finally reaches the following expressions for the crossing 
(or ordering) time. 



T 



InN 



for 6 < 1/3 
for b = 1/3 



(3i>o-l) 
(In N) 2 
3a 2 

A ~^ — for 6 > 1/3. 



(B4) 



or in terms of the original parameters 
lnN for b Q <b Q 



T 



70-1 

3(lnJV) 2 f , , 

a -i for t> = 00,, 

a(-yp-i) 

A - 2 for b a > b ,, 



(A34) 
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